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A new type of intermittent behavior is described to occur near the boundary of phase synchro- 
nization regime of coupled chaotic oscillators. This mechanism, called ring intermittency, arises 
for sufficiently high initial mismatches in the frequencies of the two coupled systems. The laws for 
both the distribution and the mean length of the laminar phases versus the coupling strength are 
analytically deduced. A very good agreement between the theoretical results and the numerically 
calculated data is shown. We discuss how this mechanism is expected to take place in other relevant 
physical circumstances. 
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Intermittent behavior is an ubiquitous phenomenon in 
nonlinear science. Its arousal and main statistical prop- 
erties have been studied and characterized already since 
long time ago, and different types of intermittency have 
been classified as types I— III [l], 0] or on-off intermit- 
tency 0, 13] . One of the most general and interesting 
manifestations of the intermittent behavior can be ob- 
served near of the boundary of chaotic synchronization 
regimes. Indeed, close to the threshold parameter val- 
ues for which the coupled systems show synchronized dy- 
namics, it is observed that the de-synchronization mecha- 
nism involves persistent intermittent time intervals dur- 
ing which the synchronized oscillations are interrupted 
by the non-synchronous behavior. These pre-transitional 
intermittencies have been described in details for the case 
of lag synchronization [B], 0, 0] and for generalized syn- 
chronization Q , and their main statistical properties (fol- 
lowing those of the on-off intermittency) have been shown 
to be common to other relevant physical processes. 

As far as intermittency phenomena near the phase syn- 
chronization onset are concerned, two types of intermit- 
tent behavior have been observed so far 
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namely the type-I intermittency and the super-l one , 
inar behavior (so called "eyelet intermittency" [13|). 

In this Letter we report that a new type of intermit- 
tent behavior is observed near the phase synchronization 
boundary of two unidirectionally coupled chaotic oscilla- 
tors, when the natural frequencies of the two oscillators 
are sufficiently different from one another. The system 
under study is represented by a pair of unidirectionally 
coupled Rossler systems, whose equations read as 

X d = -UJdUd - Z dl X r = -UJrUr ~ Z r + e{x d ~ X r ), 

y d = u; d x d + ay d , y r = uj r x r + ay r , 
z d =p + z d (x d ~c), z r =p + z r (x r - c), 

(1) 

where (x d ,y d , z d ) [(x r ,y r , z r )\ are the cartesian coordi- 
nates of the drive (the response) oscillator, dots stand for 
temporal derivatives, and e is a parameter ruling the cou- 
pling strength. The other control parameters of Eq. ([I]) 



have been set to a = 0.15 
with previous studies 
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= 0.2, c = 10.0, in analogy 
The Otv-parameter (rep- 
resenting the natural frequency of the response system) 
has been selected to be ui r = 0.95; the analogous param- 
eter for the drive system has been fixed to cu d = 1.0. For 
such a choice of parameter values, both chaotic attrac- 
tors of the drive and response systems are, at zero cou- 
pling strength, phase coherent. Furthermore, the bound- 
ary of the phase synchronization regime occur around 
e c « 0.124. 

The instantaneous phase of the chaotic signals tp(t) 
can be therefore introduced in the traditional way, as the 
rotation angle ip djr — arctan(?/d, r /x d , r ) on the projection 
plane (x, y) of each system. The presence of the phase 
synchronization regime can be detected by means of mon- 
itoring the time evolution of the instantaneous phase dif- 
ference, that has to obey the phase locking condition [l6| 



\A<p(t)\ = \<fd(t) - <Pr(t)\ < const. 



(2) 



Below the boundary of the phase synchronization regime, 
the dynamics of the phase difference A(p(t) features time 
intervals of phase synchronized motion (laminar phases) 
persistently and intermittently interrupted by sudden 
phase slips (turbulent phases) during which the value of 
|Ay>(t)| jumps up by 2ir. 

By analyzing the statistics of the laminar phases, it 
is found that the intermittent type behavior described 
in Refs. 0, [13, U, 12, 13 1 takes place only for small 
differences in the natural frequencies of the drive and 
response systems. In particular, the eyelet intermittent 
phenomenon occurs in the range u d = 0.90 -f- 0.98. As 
far as large differences in the natural frequencies of the 
drive and response systems are concerned {bj d < 0.90 and 
Lo d > 0.98), a novel type of intermittent behavior emerges 
which differs remarkably from the ones known so far. In 
the following we will describe the properties of this new 
type of behavior, that we called ring intermittency, due 
to the specific dynamical mechanism that produces it. 

Let us start with discussing the mechanism at the ori- 
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FIG. 1: (Color online) The phase trajectory of the response 
system on the (x' , y')-plane rotating around the origin when 
the coupling parameter strength is selected as (a) e = 0.126 
— the phase synchronization regime, (6) e = 0.115 — the 
ring intermittency and (c) e = 0.109 — the asynchronous 
dynamics, (d) The dependence of the phase difference A<p(t) 
on the time t for the coupling strength values e used in Fig. (a- 
c). 



gin of the arousal of the ring intermittency, and that rules 
out the scaling laws characterizing this phenomenon. It is 
well-known that there are two different scenarios for syn- 
chronization destruction in a periodic oscillator driven 
by an external force, corresponding respectively to small 
and large detunings with the external signal frequency 
(see, e.g., tutorial [l3|)- Under certain conditions (i.e., 
for the periodically forced weakly nonlinear isochronous 
oscillator), the complex amplitude method may be used 
to find the solution describing the oscillator behavior 
in the form u(t) = Re a(t)e luJt . For the complex am- 
plitude a(t) one obtains averaged (truncated) equations 
a = —iva + a — \a\ 2 a — ik, where v is the frequency mis- 
match, and k is the (renormalized) amplitude of the ex- 
ternal force. For the small v and large k the stable solu- 
tion a(t) = Ae 1 ^ = const corresponds to the synchronous 
regime, with the synchronization destruction correspond- 
ing to the the saddle-node bifurcation on the plane of the 
complex amplitude. For large frequency mismatches with 
the decrease of fc-value the fixed point (stable node) on 
the complex amplitude plane becomes sequentially a sta- 
ble focus and an unstable focus (via the Andronov-Hopf 
bifurcation). In this case the phase synchronization de- 
struction is connected with the limit cycle location on the 
complex amplitude plane (see 17f for detail). When the 
limit cycle starts enveloping the origin, the synchroniza- 
tion regime begins to destroy. Obviously, if one consid- 
ers the behavior of the synchronized periodic oscillator on 
the plane (x', y') rotating with the frequency of the exter- 



nal signal around the origin he observes the stable node 
for the small values of the frequency detuning and a cy- 
cle for the large ones, respectively. These considerations 
on the rotating plane may be made apparent by using 
the coordinate transformation x' = x r cos (fd + y r sin ip^, 
y' = — x r sm(pd + y r cos ifd, where tfd = <Pd(t) is the in- 
stantaneous phase of the drive system. 

The same effects may be also observed for chaotic os- 
cillators. Indeed, in Fig. [TJ a the behavior of the syn- 
chronized response oscillator JT]) is shown on the plane 
(x 1 , y') rotating around the origin in accordance with the 
phase tpd(t) of the drive system when the control pa- 
rameters uid and ijj r are detuned sufficiently. One can 
see that the phase trajectory on this plane looks like 
a ring. This effect arises insofar as the Rossler system 
may be considered as a noise smeared periodic oscillator 
(see, e.g., [Hj]). Therefore one observes the ring consist- 
ing of the phase trajectories, instead of the limit cycle 
that would occur in the periodic case. It is the case of 
the large control parameter mismatch that is accompa- 
nied by the ring intermittency behavior, while for the 
small parameter detuning the intermittent type-I as well 
as the eyelet intermittency are revealed. When the cou- 
pling strength e gets below the critical value e c the phase 
trajectory on the (x 1 , ?/)-plane starts enveloping the ori- 
gin (see Fig. [H b, the origin is the point of intersection of 
the dashed lines), and the phase synchronization regime 
begins to destroy, as a phase slip is observed all the 
times that the phase trajectory envelops the origin of 
that plane. As the coupling strength decreases further, 
the phase trajectory envelops the origin more often, and 
the phase slips occur more frequently. Finally, when the 
coupling strength e becomes less than et ~ 0.1097 the 
origin is inside the ring (see Fig.[TJ c), therefore every ro- 
tation of phase trajectory causes a phase slip. So, vary- 
ing the coupling strength e one observes: (i) the phase 
synchronization regime for e > e c , (ii) the intermittent 
behavior for e t < e < e c and (iii) the asynchronous dy- 
namics for e < £t when the phase slips follow each other 
at approximately equal time intervals T, the averaged pe- 
riod of the phase trajectory rotation on the (a/, y')-plane 
(see Fig. Hid). 

Let the probability that the phase trajectory on the ro- 
tating (x', y')-plane envelops the origin be p — p(e). Ob- 
viously, p = if £ > e c , p = 1 if e < s t and < p < 1 
when e t < s < £ c . 

In the case of the intermittent behavior (i.e., 
£ t < £ < £ c ) the probability of the laminar phase with 
length T to be observed is P(T) = p 2 . This period is 
determined by the difference of the main frequencies of 
the drive (fd) and response (f r ) systems and may be cal- 
culated as T l/|/r — fd\- For the control parameter set 
mentioned above, we have T w 80. It is clear that the 
probability of laminar phases with length nT to arise is 
P(nT) = (1— p)" _1 p 2 . So, the distribution of the laminar 
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FIG. 2: (Color online) The laminar phase length distribu- 
tions for the different values of the coupling strength e and 
exponential laws (j4j) corresponding to them. The theoretical 
curves are shown by lines, the data calculated for two cou- 
pled Rossler systems {1} are shown by points. The coupling 
strength values e, the probability p and the parameter of the 
exponential law k have been selected as follows: (1) e = 0.123 
(o), p = 0.06, k = -7.7 x 10~ 4 ; (2) e = 0.12 (♦), p = 0.31, 
k = -4.6xl0" 3 ; (3) e = 0.115 (■),p = 0.65, k = -OxlGT 2 ; 
(4) e = 0.11 (A), p = 0.96, k = -4.0 x 10~ 2 . 



phases with generic length r should scale as 

N( T )~p 2 (l- P r/ T -\ t > T. (3) 
Equation ([3]) may be also rewritten in the form 
N(t) = Aexp(fcr), 



t > T 



(4) 



where k — (1/T) ln(l — p), A being a normalizing coef- 
ficient. Thus, the laminar phase distribution in the ring 
intermittency obeys an exponential law, with the param- 
eter k being negative due to < p < 1. 

Let us compare the obtained theoretical law (f4]) with 
the results of the numerical calculations of the intermit- 
tent behavior of two coupled Rossler systems ([1]). In 
Fig. [2] the distribution of the laminar phase lengths are 
shown for the different values of the coupling strength 
e. In the same Figure we report also the exponential fits 
obeying the law (j4]). The value of the probability p(e) to 
calculate the coefficient k may be estimated as follows. If 
the length of the time realization under consideration is 
L (in our calculations the length of the analyzed time se- 
ries was L = 2x 10 6 time units) and N phase slips have 
been observed during this period, then the probability 
p(e) that the phase trajectory on the rotating (x',y')- 
plane envelops the origin is p = NT/L. E.g., when the 
coupling strength has been fixed as e = 0.12 (see curve 2 
and points ♦ in Fig. [2] for the distribution of the laminar 
phase lengths) we have observed N = 7, 772 phase slips 
during the analyzed time series, and therefore the prob- 
ability p(e) was estimated as p = 0.31. From Fig. [2] one 
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FIG. 3: (Color online) The dependence of the probability p 
that the phase trajectory on the rotating (x' , j/')-plane en- 
velops the origin on the deviation of the coupling strength 
from the critical value (e c — e). The points obtained numeri- 
cally has been shown by symbols 0, the linear approximation 
p = a(e c — e) (where a = 70) is shown by a solid line. The crit- 
ical values s c and e t of the coupling parameter are shown by 
arrows. Note, equation (0 is applicable only in the coupling 
strength range et < e < e c . 



can see an excellent agreement of the numerical data with 
the theoretical law ((4]) for the whole range of coupling 
strength values St < s < s c where the ring intermittency 
takes place. 

Let us now derive the dependence of the mean length 
(t) of the laminar phases (i.e., the averaged time inter- 
val between two successive phase slips) on the coupling 
strength e. From equation (U|) one can easily obtain 
the relationship between the mean length (r(e)) and the 
probability p = p(e) 



(r)=T 



T 



T 



ln(l-p) 



(5) 



We have numerically observed that, in the coupling 
strength range e t < e < e c where the ring intermittent 
behavior is observed, the value of the probability p is 
directly proportional to the deviation of the coupling 
strength e from the critical value e c , i.e., 



p(e) - (e c - e) . 



(6) 



Indeed, in Fig. [3] the dependence of the probability p on 
the deviation of the coupling strength from the critical 
value (e c — e) is shown. The probability p(e) for each 
value of e has been calculated in the same way as it was 
described above. 
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FIG. 4: (Color online) The mean length (r) of the lami- 
nar phases vs. e calculated numerically (0) for two coupled 
Rossler systems {l}, the theoretical curve Q (red solid line 
1, applicable only in the range e t < e < e c shown by arrows), 
and the asymptotic (e < e t ) value of ([5]) (red dashed line). 
The analogous curves (filled dots, blue solid line, blue dashed 
line) refers to the case of two coupled generators with tun- 
nel diodes of Ref. In this latter case, arrows point to 
£t2 = 0.083 and e C 2 = 0.111, delimiting the validity range of 
Eq. ©. 



Since the probability p on the coupling strength in 
the range e t < e < e c relates to the coupling strength as 
p(e) = (e c — s)/(e c — Et), one easily obtains that the de- 
pendence of the mean laminar phase length on e has to 
scale in the form 



<T(e)> = T 1 - In 



the theoretical equation 0. Again, one can see a per- 
fect agreement of the theoretical curve with the calcu- 
lated points. Below the threshold e t , the calculated mean 
length of the laminar phases is compared with the asymp- 
totic value (r) = T (red dashed line). In order to show 
that our analysis is not limited to the chaotic system fU , 
we report in the same Figure the analogous curves (filled 
dots, blue solid line, blue dashed line) obtained for the 
case of two unidirectionally coupled generators with tun- 
nel diodes described in Ref. [l9| , where Eq. (O is valid in 
the range e t i < £ < e c i (e t 2 — 0.083 and e c2 = 0.111). Fi- 
nally, we would like to stress that the same intermittent 
scenario can be observed in system (1) for fixed (large 
enough) values of coupling strength, when varying the 
parameter mismatch. 

In conclusion, we have reported for the first time a 
new type of intermittency behavior occurring at the on- 
set of phase synchronization regimes of two unidirection- 
ally coupled chaotic oscillators with sufficiently detuned 
natural frequencies. Such a type of ring intermittency 
differs remarkably from all the other types of intermit- 
tency known so far. It may be observed in a certain range 
of coupling parameter strengths, where the distribution 
of the laminar phase lengths obeys an exponential law. 
The theoretical equation for the dependence of the mean 
length of the laminar phases on the coupling strength 
has also been given, and are in perfect agreement with 
the numerically obtained data. Though the character- 
ization of the new intermittent process has been here 
explicitly derived at the boundary of phase synchroniza- 
tion of chaotic systems, we expect that the very same 
mechanism can be observed in many other relevant cir- 
(7) cumstances, as e.g. laser systems [11], or in the case of 



Notice that Eq. ([7]) describes scaling properties for 
the laminar periods during the ring intermittency that 
are completely different from those typical of type-I inter- 
mittency characterizing the transition to phase locking of 
periodic oscillators in the presence of noise [l6[ , as well as 
from those arising from the super-long laminar behavior 
(or eyelet intermittency) characterizing pre-transitional 
stages of slightly mismatched chaotic oscillators 11 , 13[ . 
Obviously, equation ([7]) is correct only in the coupling 
strength range e t < e < e c , whereas §5$ may be used both 
below Et and above s c . The mean length (r) of the lam- 
inar phases is about T for values of e is below the crit- 
ical value £t (the asynchronous regime where the phase 
slips follow one another at approximately equal time in- 
tervals), and is infinity in the case e > e c . Note also, 
that lim e ^ +et (r(e)) = T and lim £ _ > _ Ec (r(e)) = +oo in 
perfect agreement with the system's behavior at the two 
boundaries of the intermittency phenomenon (e t ;e c ). 

Finally, in Fig. [4] the dependence of the mean lami- 
nar phase length (r) on the coupling strength e is shown. 
The points (0) correspond to the calculated mean lengths 
of the laminar phases and the red solid line 1 reports 



the interaction between the main rhythmic processes in 
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the human cardiovascular system 
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